#include "Possion.h"

int main(int argc, char* argv[]){
  /*
  Point<double,2> P({1,2},0,0);
  std::cout << "exact solution:" << function_u(P) << std::endl;
  std::cout << "action:" << function_f(P) << std::endl;
  */
  EasyMesh Mesh("data/D");
  VectorXd Solution = Possion_Easymesh(Mesh,1,1e-15);
  //double useless_err = useless_error(Mesh,Solution);
  //std::cout << useless_err << std::endl;
  double prior_err = prior_error(Mesh,Solution,3);
  std::cout << prior_err << std::endl;
  //plots_matlab(Mesh,Solution);
  return 0;
}
